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The photon detection efficiencies of multi-pixel Geiger-mode avalanche photodiodes 
manufactured by different producers are estimated. A new fit method of the response 
spectra to low-intensity light, taking into account after-pulse and cross-talk effects is 
proposed to yield the initial number of photons. The value of photon detection efficiency 
is calculated using a calibrated photodetector as a reference. 

1 Introduction 

Multi-pixel Geigcr-modc avalanche photodiodes (G-APDs) arc solid-state photodctcctors 
based on a technology developed since early 1990's [IJ[2l[3l[4j[5j[6l[7l[8l[9]. Their gain 
is typically 10 6 per photoelectron. Their main features are: the small size, the efficiency 
comparable to the one of vacuum photomultipliers, the insensitivity to magnetic field, a low 
bias voltage and a reasonable price. Currently, the G-APDs find an application in many 
fields of physics, High-Energy, Neutrino, Astroparticle and Medicine. 

This paper describes the results of a comparative study performed in the same experi- 
mental conditions, on different G-APDs, to allow an optimal choice with a view to a given 
application. We refer to elsewhere for a complete description of measurements and analysis 
methods [10] . 

2 General principle of the G-APD operation and its characteristics 

A G-APD consists of a large number of identical microcclls (silicon diodes, pixels) with a 
common anode. The microcclls arc located on a common substrate with a typical size of 
~ 1 x 1 mm 2 . Once photoelectron are produced an electron avalanche multiplication may 
occur under a reverse bias, Vbi as , above the breakdown voltage, Vbrd and the charge induced 
during the avalanche discharge is proportional to the overvoltage (V OVC r = Vbias -Vbrd)- The 
number of fired pixels is proportional to the initial number of photoelectrons, as long as 
this is on average less than one per pixel. This implies a low-intensity light source in the 
measurement setup. The gain (g) is determined by the charge released in one pixel. 

Other characteristics of the G-APD response to light deserve a few words. 

After-Pulses (AP) appear when the quenching (mechanism to stop the avalanche process 
by temporally reducing the cell voltage) doesn't drain all the charge in the sensitive area, 
and the cell fires again shortly after the original pulse. 

Cross-talk takes place when in a Geiger discharge, some of the electrons generated in 
the avalanche process reach the adjacent pixel, where another Geiger discharge might be 
triggered. This results in a pulse with a larger amplitude (nearly a factor two) and in a 
distortion of the linear response of the device. 
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Every detector can only convert a certain percentage of incident photons in signals. This 
G-APD efficiency is given as product of three factors: 

PDE = QE X e gcom X e G eiger, (1) 

where QE is the G-APD quantum efficiency, function of the incident photon wavelength, A w ; 
fgoom is a geometrical factor (fill factor) indicating which fraction of the device is sensitive 
to photons; eooiger is the probability to trigger a Geiger discharge and is function of Vbias- 
PDE decreases with the pixel size. 

Then, the G-APD response can be written as a product of few factors: 



A = N! y n xPDExgx (1 + e) x (1 + Pap), (2) 

where N™ is the number of incident photons; g, the gain; e, the cross-talk probability and 
Pap, the after-pulse probability. It should be noticed that in the above equation, g, PDE, 
e and Pap are all increasing with Vbias, implying a complex dependence of the G-APD 
response on the bias voltage. 



3 Measurement setup 



Figure [T] shows the scheme of the PDE measurements. The light from a Light-Emitting 
Diode (LED) operating in a pulse-mode with a Waveform Generator (AGILENT 33220A) 
was delivered to an optical filter. The latter corresponded to the peak wavelength of the 
LED and had a bandwidth of ± 3 nm. By changing the LED, a few wavelengths in the 380 
-T- 650 nm range were investigated. For low intensity sources (380 nm, 400 nm, 565 nm) the 
LED light was not filtered. 

The filtered light was routed to a light-tight thermostabilized box with inside two photo- 
detectors: the G-APD and the reference photodetector. The box temperature was controlled 
at ± 0.1 °C during the data taking. As a reference photodetector a PhotoMultiplicr Tube 
(PMT) H5783P produced and calibrated (for this reference) by HAMAMATSU was used. 
For each wavelength A w from the efficiency of the reference photodetector for the peak 
wavelength, as provided by the manufacturer, is derived the efficiency of the reference device, 
e evai' f° r our u §ht source, taking in account the filter bandwidth or the LED full width half 
maximum. 

The signal from G-APD is read out with 
a charge-sensitive preamplifier and digitized 
with an integrating ADC (CAEN QADC 
V792N). The LED pulse of about 6 ns du- 
ration and the ADC gate of about 65 ns 
width were synchronized by means of a com- 
mon trigger. The integration time was long 
enough to contain most of G-APD signal. In 
case of the reference photodetector the di- 
rect signal from PMT was amplified with an 
external amplifier (NIM ORTEC 450) and 
digitized with the same ADC. Due to a high 
sensitivity to voltage variations, the G-APD 
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Figure 1: General scheme of the 
surcmcnts. 
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bias voltage was provided by a high precise 

and stable in time voltage source (Keithley 6487/E). 



4 Samples 

The measurements refer to the samples listed in Tab. [TJ The following operating conditions 
were chosen: 

• HAMAMATSU produced Multi-Pixel Photon Counter S10362-11-025U [H] operated 
at Vbias = 71.32 V corresponding to V ovcr ~ 2.3 V. 

• CPTA produced SiliconPhotoMultiplier operated at Vbias = 32.5 and 32.7 V corre- 
sponding to Vovcr — 2.3 V and 2.5 V. The parameter dependence on voltage has been 
studied in the range between 31.5 and 33.3 V. 

• IRST produced SiliconPhotoMultiplier operated at Vbias = 32 V corresponding to 



The temperature at data taking was T ~ 23.7 °C stabilized at ± 0.1 °C during each mea- 
surement and with a spread ± 0.3 °C during all period of data taking. 



Sample 


Type 


Photosensitive area 
[mm 2 ] 


Number of pixels 


Pixel pitch 
[/zm 2 ] 


HAMAMATSU 


S10362-11-025U 


lxl 


1600 


25x25 


CPTA[1] 


143 


1.028[2] 


556 


43x43 


IRST [3] 


2007 prod 


lxl 


400 


50x50 



Table 1: Specification of measured samples. From left: Manufacturer, type, photosensitive 
area, number of pixels and their pitch. [1] Sample kindly provided by Prof. M. Danilov, 
[2] Sensitive area has octagonal shape, [3] Sample kindly provided by Dr. C. Picmontc 
(Fondazionc Bruno Keller). 



5 Measurement strategy and fit procedure 

As described in Sec. [3] the PDE measurement strategy is based on the comparison of the 
effective number of photons detected from reference detector and G-APD. 

To accomplish this task, noise and light (signal) response spectra of G-APD and PMT 
to a low number of photons have been measured and fitted, as described in the following. 

The number of photons, n, detected per light pulse is expected to be Poisson distributed, 
with mean value A. This distribution has to be convoluted with a gaussian distribution 
describing the experimental resolution. The charge x measured by the ADC is expected to 
have a frequency distribution, N(x): 



n>0 



N(x) = N x J Poisson(A) ® Gauss(^ n , cr n ) j , 
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or more explicitly: 



N(x) = N x e _A 



A" 



1 , -(x- Mn ) 2 

— x — exp( — 

n! cr n 2o-,l 



(•3) 



In the ideal case n photons correspond to a charge fx n = fiQ + n x g, where fio is the 
pedestal position and g is the gain factor. In presence of a finite resolution the distribution 
of charge follows a gaussian distribution around each /i n with a width o~ n , o~\ = (7q + n x o~\ , 
where (To is the electronic noise (pedestal) width and o\ is the single photon width. N is a 
normalization factor. 
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Figure 2: Reference PMT: noise (left) and signal (right) spectra. The signal is obtained by 
short, low-intensity LED flashes at A w = 450 ± 3 nm. The continuous lines are the fit results 
(see text). From the fit the mean number of detected photons is obtained, N™ T = 1.71±0.02 
( e evli= 20 - 33 tols^, see text). Binwidth^l. 



Reference PMT fit procedure In the PMT case, it should be noticed that due to additional 
electrical noise caused by the external amplifier, the no-light spectrum (pedestal) is described 
as superposition of three gaussian peaks, two of which arc considered to be symmetric. 
Figure [2] (left) shows the pedestal fitted with the sum of three gaussians. Figure [2] (right) 
shows the signal fitted with Eq. [3] modified to include the noise shape. 

Using the mean number of detected photons, A = N™ T , as obtained from the above fit, 
and the PMT efficiency, the mean number of photons delivered by the optical system to the 
photodetector surface per one LED pulse was estimated. 

In the example of Fig. [U at A w = 450 nm, it is N™ T = 1.71 ± 0.02. Taking into account 
the efficiency for this A w , e^i = 20.33 lois %, the mean number of photons impacting 
PMT is obtained. This is the first step for the PDE calculation. 

G-APD fit procedure Figure [3] shows typical G- APD HAMAMATSU response spectra, on 
the left for no light (noise) and on the right for low- intensity light (signal). 

The noise spectrum shows beyond the first peak (pedestal) other peaks which are due to 
thermogeneration. Since their contribution to the total noise is about 2% of the pedestal, 
the thermogeneration effect can be neglected in the signal fit. Moreover since the probability 
of two or more thermogeneration counts within the gate is negligible, the cross-talk effect 
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Figure 3: HAMAMATSU G-APD: noise (left) and signal (right) spectra. The signal is 
obtained by short, low- intensity LED flashes at A w = 450 ± 3 nm. Each peak corresponds 
to a number of detected photons. Binwidth=l. 



can be estimated from the ratio of the numbers of events in the 3 and 2 nd peak. This 
value (~ 20%) depends on the gate length (in our case 65 ns) and on the operation voltage, 
Vbias (Sec. E). 

In the signal spectrum the first peak (pedestal) corresponds to the noise, the second 
one is the G-APD response when exactly one photon is detected (one pixel is fired). It is 
impossible to identify which one of the pixels was fired, since they are all connected to the 
same output and they have similar responses. If n pixels are fired, the sum of n charges is 
recorded at the position of the (n + l) th peak, fx n = + n x g. 

These peaks, as shown in Fig. [3] (right), are equidistant, at a distance determined by 
the gain factor, g. The enhancements which show up between the peaks, were attributed to 
after-pulses (AP) and were considered in the fit. 

To identify the most appropriate function describing this effect the main peaks, each 
fitted to a gaussian distribution, were subtracted from the spectrum of Fig. [3] (right). As a 
result, a sequence of peaks, Fig. |4] (left), was obtained, each one well described as a sum of 
two gaussians (for this G-APD and for the chosen gate length). This suggested to write a 
new fitting function P(x) as follows. 

If P„ is the probability to get initially n cells fired (the poissonian n term in Eq. [3|), 
the probability to record an x charge, P(x)=N(x)/N can be written as sum of three terms 
representing respectively the probability of an x charge when n cells are fired 1) without 
AP, P£ oAP ; 2) when only one AP occurred, P„ P1 ; 3) when two APs occurred, P„ P2 : 

P(x) = P n° AP x Gauss(ju„, cr„) + P AP1 x Gauss(^„ + 5i,ai) + P„ P2 x Gauss(/z n + S 2 , a 2 ), 

n 

(4) 

where are the distances of the first and second after-pulse gaussian from the nearest 
main peak and 01,2 the widths. The probabilities PJJ oAP , P„ P1 and P„ P2 are all functions 
of a single parameter Pap, the probability to get an AP from one cell, and can be written 
as: 
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Figure 4: Residual distribution after subtraction of best-fitted gaussians from the signal 
spectrum (left). The signal spectrum fitted with the after-pulse corrections (right). From 
this fit, the probability, P°, to get n cells fired and the probability, Pap, to get an after-pulse 
from one cell are derived (in this example Pap = 0.16 ± 0.01). Binwidth=l. 



P7 AP = p° x (i - p AP ) n , 

P AP1 = P° x (1 - (1 - P AP ) n ) x (1 - P AP ) n , 
P AP2 - P° x (1 - (1 - P AP ) n ) x (1 - (1 - P AP ) n ). 

If the probability of after-pulses is zero, the Eq. |4]is equivalent to Eq. GS 

Equation |4] was fitted to the spectrum of Fig. [3] (right) and the fit result is shown in 
Fig. H] (right) as a continuous line. Among the free parameters, we obtain the P° values as 
displayed in Fig. [5] (black dots) and the after-pulse probability, Pap = 16%. 

The P° values should be distributed according to the Poisson statistics with a mean 
number of photons detected, N^ _APD (corresponding to the Poisson's distribution parameter 
A). However this distribution is distorted by the presence of cross-talk. If the cross-talk 
probability is e ^ 0, then the n=l bin (2 nd peak in Fig. |4] (right)), contains only events with 
one initially fired cell without any cross-talk, with a probability P(l)= P?(l — £)■ The n=2 
bin (3 rd peak) is filled either when two cells are initially fired without cross-talk or when one 
cell is fired accompanied by a second cell fired by cross-talk: P(2) = P°(l — e) 2 + P?e(l — e). 
Finally the probability to observe n fired cells P(n), can be written: 

n / 1 

-n „ : /n — 1 



P(n)=^P J °(l-ey £ "-j . (5) 

j=i VJ ' 

The data in Fig. [5] (left) were fitted using Eq. [5j the red line represents the result of the fit. 
For these data, the value obtained for the cross-talk probability was e = 0.20 ± 0.01; the 
value of the mean number of detected photons is N^ _APD = 1.55 ± 0.02. In general, for any 
pair of measurements on G-APD sample, and on the reference PMT, we have: 

|^G-APD 

PDE = N PMT . e re f ( 6 ) 
iN 7 fc eval 
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Figure 5: Right: distribution of probability, P, to observe N 7 photons. The data, black dots, 
are the P° obtained by fitting the signal spectrum (Fig. [4] right). The red line represents the 
fit of these data with Eq. [5] In this example, the fitted values are: the mean number of G- 
APD detected photons, N^ APD = 1.55±0.02, and the cross-talk probability, e = 0.20±0.01. 
Binwidth=l. Left: The photon detection efficiency, PDE, for the HAMAMATSU G-APD 
as a function of the light wavelength, . 



With the mean number of photons detected by G-APD, N^" APD and by PMT, the 
reference detector, N PMT , it is possible to evaluate the photon detection efficiency. In our 
example, the PDE value is 18.36 ± 0.26 % (statistical error only). 

The final value of PDE is obtained as mean of various measurements at different light 
intensities, with different fiber configurations, taking into account the FWHM of the light 
source and the filter bandwidth. 

The procedure and results described in this section for G-APD HAMAMATSU were as 
well applied to the others samples, with minor modifications. In the following sections these 
results arc discussed. 



6 Photon detection efficiency: Hamamatsu sample 

Based on the fitting procedure of Sec. [5J it is possible to calculate for each measurement the 
corresponding PDE value. The measurements have been made at different light intensities 
and with different optical layouts to evaluate possible sources of systematics. The systematic 
error due to optical contact, excluding a statistical contribution, is ~ 1% in absolute value. 

At this point, it is possible to evaluate the photon detection efficiency of the device as the 
mean of all available measurements. The results are shown in Fig. [3] (right). The horizontal 
error bar is originated from the light source FWHM or filter bandwith; its impact on PDE 
evaluation gives the vertical error bar. 
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Figure 6: The noise (left) and signal (right) spectra obtained with the CPTA G-APD illu- 
minated by short, low intensity LED flashes at X w = 600 ± 3 nm. Each peak corresponds 
to a number of detected photons. Binwidth=l. 

7 Photon detection efficiency: CPTA sample 

The noise and the signal spectra obtained with the CPTA sample, as shown in Fig.[S]for X w 
= 600 ± 3 nm, reflect the difference on the construction characteristics. With respect to the 
HAMAMATSU device Fig. [6] (left) shows that with the chosen gate the thermogeneration is 
at a level of ~ 30% , and Fig. [6] (right) shows that the after-pulses have lower impact on the 
single photon spectra. The cross-talk effect is much lower. The fit procedure adopted for 
the CPTA device is slightly different from that described in Sec. [SJ The thermogeneration 
contribution was not included in the fit, the AP terms ( Eq. 3} were not directly related to 
an after-pulse effect but accounted also for the termogeneration and other effects. 

The PDE results for other wavelengths are shown in Fig. [7] (left). The horizontal error 
bar is originated from the light source FWHM or filter bandwidth, the vertical one from 
its impact on the PDE evaluation. The measurements have been recorded at Vbi as =32.5 V 
for X w = 450 nm, 565 nm, 600 nm and Vbias=32.7 V for X w = 500 nm and 600 nm. The 
measurement at X w = 600 nm has been recorded at both voltages. 

7.1 Parameter dependence on overvoltage: PDE 

The photon detection efficiency, PDE ro i, relative to the one at the reference voltage as 
a function of Vbias has been recorded at different wavelengths. From a straight line fits 
(PDE rc i = po + Pi Vbias) these distributions the values for the angular coefficients, pi[V _1 ] 
are derived and the PDE value as a function of X w at various values of Vbias is obtained, 
Fig. [7] (right) . The validity of the method is tested with the measurements at 600 nm, taken 
at two different bias voltages. 

7.2 An improved fit procedure 

The CPTA G-APD measurements described above have two peculiar aspects that deserve 
more attention: the thermogeneration effect is not negligible (Fig. |6] (left) and the after-pulse 
effect requests a dedicated treatment. We refere to [10], for a detailed description of the 
method. 
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Figure 7: Left: the photon detection efficiency, PDE, for CPTA G-APD as a function of the 
light "wavelength, A^, . The measurements at 450 nm, 565 nm, 600 nm have been recorded at 
Vbias=32.5 V, those at 500 nm, 600 nm at Vbi as = 32.7 V. The X w = 600 nm measurements 
at the two bias voltages are indistinguishable. Right: the photon detection efficiency, PDE, 
for the CPTA detector as a function of the wavelength, A™, as measured (red points) and 
extrapolated (green, blue, magenta points). 



As a first step the termogeneration spectrum shape is determined as a deviation of the 
noise spectrum, Fig. [5] (left), from a gaussian. A combination of two exponential functions 
results into an appropriate description of this effect. 

As second step the fit of the signal pedestal is not performed with a simple gaussian 
as above, but includes also the fit function determined in the first step and describing the 
termogeneration. 

Further improvement in the signal fit is 
achieved, if the after-pulse effect is modeled 
for a single fired pixel, as suggested in Sec. 
[21 The results of this improved method are 
shown in Fig. [8] 

Using this fit procedure it is possible to de- 
rive a value of the cross-talk probability and 
the average obtained on measurements per- 
formed at many light intensities is < 1.7%. 

Similar measurements were recorded 
with a IRST device. Again, the noise and 
signal spectra reflect the difference on con- 
struction and design. From these spectra 
we have concluded that the specimen used 
wasn't suitable for PDE measurements. 
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Figure 8: The CPTA G-APD signal spectrum 
(at A w =450 ± 3 nm) fitted with the improved 
fit procedure, Sec. 17.21 
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8 Conclusion 



Measurements of G-APD response to low intensity light were recorded to determine the main 
characteristics of three different samples, HAMAMATSU (S10362-11-025U), CPTA, IRST. 
An accurate method to fit the response has been realized, to evaluate photon detection 
efficiency, cross-talk, gain, and after-pulse probability. A calibrated reference detector has 
been used for the photon detection efficiency measurements in a wavelength range between 
380 and 650 nm, see Ref. [TO]. 

The measurements were carried on in the same experimental condition (setup, thermal 
insulation...) to compare the properties of various samples. This information is an important 
issue in the near future to identify the most suitable device for a particular application in 
high-energy physics calorimeters, astroparticle physics and medical environment. 

The PDE measurements provided are consistent with those quoted from the manufac- 
turers with a caveat due to the different method of measurement (the current method, very 
often used by manufacturers, includes after-pulse and cross-talk effects). 
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